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Abstract. A dyadic shell model for the Navier-Stokes equations is studied in 
the context of turbulence. The model is an infinite nonlinearly coupled system 
of ODEs. It is proved that the unique fixed point is a global attractor, which 
converges to the global attractor of the inviscid system as viscosity goes to zero. 
This implies that the average dissipation rate for the viscous system converges to 
the anomalous dissipation rate for the inviscid system (which is positive) as vis- 
cosity goes to zero. This phenomenon is called the dissipation anomaly predicted 
by Kolmogorov's theory for the actual Navier-Stokes equations. 



We consider the Navier-Stokes equations for the motion of a three-dimentional 
incompressible viscous fluid: 



Here u denotes the velocity vector field, p the pressure, / an external force, and 
u the viscosity coefficient. The role of the nonlinear term in (11.11 ) is critically im- 
portant in the theory of turbulence where a basic principle is a cascade of energy 
from large scales, through the so called inertial scales, to very small dissipative 
scales. Transfer of energy through these scales is achieved via nonlinear interac- 
tions between the modes in the Fourier space. This subject is the topic of extensive 
study in the experimental, numerical, and analytical literature (see, for example, 
Frisch fTl ], Eyink and Sreenivasan [TU]). Important seminal work in the modern 
theory of turbulence was performed by Kolmogorov and his school in the mid 20th 
century. However, rigorous mathematical proofs of Kolmogorov's laws remain to 
be obtained. 

Kolmogorov predicted that the energy cascade mechanism in fully developed 
three-dimensional turbulence produces a striking phenomenon, namely the persis- 
tence of non-vanishing energy dissipation in the limit of vanishing viscosity. This 
behavior, called the "dissipation anomaly", can be described as follows 



1. Introduction 



(1.1) 




Vp — i^Au + /, 



V ■u = 




as 1/ — > 0, where is a solution to (11.11 ). 
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Onsager conjectured that sufficiently rough solutions to the Euler equations (i.e., 
(II. lb with = 0) can exhibit turbulent, or anomalous dissipation. More precisely, 
if the Holder exponent h of the velocity is greater than 1/3, then energy is con- 
served, however, this ceases to be true if h < 1/3. For recent results concern- 
ing Onsager's conjecture see, for example, Eyink [9], Constantin, E, and Titi [6], 
Duchon and Robert [8], Cheskidov, Constantin, Friedlander, and Shvydkoy El. 

Partially because of the difficulty of proving mathematically rigorous results in 
turbulence theory, a number of "toy" models that preserve some features of the 
nonlinearity of the fluid equations have been proposed and studied by physicists 
and mathematicians. These include the so called shell models of the energy cas- 
cade, where the nonlinearity of the 3D NSE is simplified by considering only local 
interactions between scales. In this paper we study one of the original shell models 
introduced in the context of oceanography by Desnyanskiy and Novikov Q. This 
model, referred to as the dyadic model, can be written as the following infinite 
system of coupled ordinary differential equations: 

(1.2) ^aj + 1^2^^ aj - 2'^^~^^a]_, + 2'^ajaj+i = fj, j = 0, 1, 2, . . . , 

where a_i = 0, c is a positive parameter, and ^aj represents the total energy in 
the frequencies of order 2^ . For convenience we chose the force / so that /o > 
and fj = for all j > 0. 

This model has been analytically studied by Katz and Pavlovic |[T3l . Cheskidov 
ifn . Onsager's conjecture for the inviscid model (i.e., (11.21 ) with u = 0) was proved 
in Cheskidov, Friedlander, and Pavlovic [3, 5|, where it was shown that the invis- 
cid system exhibits anomalous dissipation and the unique fixed point is a global 
attractor. 

Consider (11.11) on the whole space = and define Sju as follows: 

SjU = u*T~^{ip(-2-^)), 

where is a smooth nonnegative function supported in the ball of radius one 
centered at the origin and such that V^(^) = 1 for < 1/2, and T is the Fourier 
transform. Then we define the energy flux due to nonlinear interactions through 
the sphere of radius 2^ (see [4}) as 

Hj = — u ■ VS'^u ■ u dx. 

Using the test function S'^u in the weak formulation of the Navier-Stokes equations 
we obtain 

(1.3) ~S]u = -Uj-u\\WS]u\\2. 

Recently, Cheskidov, Constantin, Friedlander, and Shvydkoy lU obtained the fol- 
lowing new bounds on the nonlinear term in (II. lb : 



(1.4) |n,-| 



u ■ VSj u ■ udx 



< V 2"ilJ-^l2^||uiii3 
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where uj is a Littlewood-Paley piece of u defined as 

Uj = Sj-^-iu — SjU. 

This estimate employing the Littlewood-Paley decomposition provides detailed in- 
formation concerning the cascade of energy through frequency space. More pre- 
cisely, it shows that the energy flux Ilj through the sphere of radius k is controlled 
primarily by scales of order k. 

Recall that Bernstein's inequality can be stated as 

(1-5) \\nj\U<2^^/P-'/'i^^\\u,\\p, l<p<q. 

Now if we define aj = ||nj||2, then using Bernstein's inequality with p = 2 and 
g = 3, we obtain 

(1-6) a^j<hj\\l<2^^/'al 

Motivated by (11.41 ) and (11.61) . we model the flux in the following way: 

(1.7) n,- = 2'^a]aj+^, 

where the scaling parameter c G [1,5/2]. Here the bounds for the scaling pa- 
rameter c are determined by Bernstein's inequality (11.51 ). with the upper bound 
corresponding to saturation of the inequality. In a turbulent flow it is expected that 
the degree of such saturation could vary giving a rise to a phenomenon known as 
intermittency. Motivated by (11.31 ) we write 




Subtracting (11.91 ) from (11.81) gives 



(1.10) ~a] = n,_i - n,- - i^2^a], 

which from the definition of the flux results in the model system (11.21) . 
In this context, the energy E and the Sobolev norms are defined as 

^ ^ oo oo 

j=0 j=0 

In this present paper we prove the following results for the viscous model ( 11.21 ). 
In Section 2 we study a steady state a, which has a property that aj is monotonic. 
This property is proved for 3/2 < c < 5/2. Hence the rest of the results in this 
paper are valid only for this range. In Section 3 we study the long-time behavior 
of solutions to (fOl ) with c > 3/2 and initial data a(0) G aj{0) > for all 
j > 0, and prove that the fixed point a is a global attractor. Moreover, a converges 
to the fixed point of the inviscid system (which is a global attractor of the inviscid 



4 



ALEXEY CHESKIDOV AND SUSAN FRIEDLANDER 



system) as u ^ 0. This allows us to conclude that the average dissipation rate 
for the viscous system converges to the anomalous dissipation rate for the inviscid 
system (which is positive) as — s- 0. 

Acknowledgments. A.C. was partially supported by NSF grant numbers DMS 
0807827. S.F. was partially supported by NSF grant numbers DMS 0803268 and 
DMS 0503768. 



2. The hxed point 
In this section we study steady state solutions a to (11.21 ). We rescale the variables 

by 

to obtain the system of equations for steady states 

^2 1) 4-1 - ^J-^^+i = '"2*'^^ ' i = 1> 2, . . . 

-AoAi = hAq - 1, 

where fj, = i/2~^/^/q and (3 = 2(1 — c/3). For the inviscid model studied 
in 131151, there is a unique fixed point with an explicit expression, namely {Aj = 
1}. In the case > we cannot solve (12.11) explicitly. However, we study fixed 
points using the monotonicity property (Theorem I2.2I ). which we prove following 
the analysis given by Hey wood \TT\. He treated the particular case of (12.11 ) where 
(3 was set to 0. In Section 3 we will show that the fixed point is a global attractor, 
i.e., it is unique. 

We are only interested in finite energy solutions a G P, which translates to 
A € H"^/^. The proof of the existence of a solution to (12.11 ) follows from standard 
Navier-Stokes techniques in which fixed point arguments are used for truncations 
of the system. Standard arguments also show that all these solutions are in H'^ for 
all s > —5/6. Moreover, we have the following 

Lemma 2.1. Let A G H^^^^ be a solution to (12.11 ). Then Aj > Ofor all j and 
where J is such that n2l^^ = 1. 

Proof. Since A G for all s > 0, we have that Aj — > as j — > cx). Then (12.11 ) 
implies that 

oo 

^,Y,2P^ A] = A,, 

3=0 

and 

oo 

lxY,2P=A] = A]_^Aj, J>0. 
j=J 

Hence Aj > for all j. 
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Now note that (12.11 ) gives 
(2.2) ^[■^'+^ ^ < i2-/5([j]+fe) < 2-/3^ 

□ 

Note that since Aj — > as j — > cxd, Lemma [ZT] implies that Aj decays super- 
exponentially. This result does not depend on the monotonicity of the sequence 
{Aj} and hence it holds in the whole range c € [1,5/2]. 

Theorem 2.2. Every solution A € H~^/^ of (12.11) with c > 3/2 is monotonic, i.e., 
Aj^i > Ajforallj > 0. 

Proof. Let hj = Aj — Aj^i. Then (12.11 ) gives 

^2 3) ^J+i = -^J- - ^2*' - ^-4i-l/^J, j > 0, 

/ii = -ho - fj. + I/Aq. 

We prove that hj < for all j by contradiction. Assume that hj > for some J. 
Then hj+i < -^il^'^ < 0, i.e., Aj+i/Aj > 1. Then ([23]) implies that 

/ij+2 > 2;u2^'^ - ^2^(-^+i) = 1x2^^(2 - 2^). 

Since c > 3/2, we have that /? < 1. We conclude that /ij+2 > 0. Iterating this 
process we obtain 

and 

^J+2A: > 0, 

for all k>0. Then 

which contradicts the fact that H'^ norm of A are finite for all s > 0. □ 
Lemma 2.3. Let Abe a fixed point. Then 



lim 



for every j > 0. 

Proof. The equations (12.11 ) read 



(2.4) 



1 

^0 



2 

^■+1 = ^-/^2^^ j = 0,1,2,... 

We will proceed by induction. Suppose that Aj^i — > 1 as /i ^ for some j > —1. 
We will show that — 1 as ^ — > 0. Assume the contrary. Then we can pass to a 
subsequence /i„ — as n — > cx), such that either 

lim sup < 1, or liminf Aj > 1. 
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First, assume that limsup^^_^o < 1- Then (12.41 ) implies that lim inf^.^^n > 
1, which contradicts the monotonicity of A (Theorem I2.2I ). Now assume that 
liminf^^^o > 1- Then (12.41 ) implies that 

limsupylj+i < 1, and liminf > 1) 

which again contradicts the monotonicity of A. □ 

The following property of a fixed point A will be used in Section 3. 
Lemma 2.4. Let A be a fixed point. Then there exists 7 G (0,1), such that 

/or all fj, > and j > 0. 

Proof. The monotonicity of A implies 

(2-5) < t4^- 

Aj + 

From (12.11 ) and monotonicity we also obtain 

(2.6) ^2+1 - //2^(^+l)^,+2 = Aj+2Aj+3 < ^■+2- 

Hence 



(2.7) A,+2 > -i/i2'3(j'+2) + i^/x222/30-+2) + 4^2_^^. 
We define 

From ( 12.11 ) and the positivity of each Aj we conclude 
Hence, 

(2.8) yz < 2^-\ 
Substituting (1277] ) into ( |23] ) gives 



5i(M) < 



1 — Z + Y^y2 _|_ -,2 

subject to constraints (12.81) . < y < 1, and < z < 00. Hence 



y2 

y - 2/3-1 _^ y/y4 _^ 22{/3-i) 
Since || > 0, /i attains its maximum at y = 1. Thus 

< V < (1 - 7)2"^/', 

' 1 - 2/5-1 + ^1 + 22{/3-i) 

provided (3 < 1, i.e., provided c > 3/2. □ 
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3. The global attractor 

In this section we study the long-time behavior of solutions to the time depen- 
dent system We study the viscous dyadic model 

(3.1) j^a, - 2^(^-i)a2_i + 2^^a,a,+i + v2^aj = j = 0, 1, 2 . . . , 

where a_i = 0. Here c € (3/2, 5/2], and i/ > is the viscosity. The initial data is 
assumed to be a(0) € aj(0) > for all j. 

Definition 3.1. A solution of ( 13.11 ) on [T, oo) (or (-oo, oo), if T = -oo) of (13.11 ) 
is an /^-valued function a{t) defined for t G [T, oo), such that aj S C^([T, oo)) 
and aj{t) satisfies ( 13.11 ) for all j. 

Note that if a{t) is a solution on [T, oo), then automatically aj G C°°([T, oo)). 
The following theorems were proved in lH. 

Theorem 3.2. For every G with > there exists a solution of (13.11 ) with 
a(0) = a°. Moreove?; aj{t) > Ofor all t > 0. 

Theorem 3.3. Let a{t) be a solution to (13.11 ) with aj(0) > 0. Then a{t) satisfies 
the energy inequality 

(3.2) \a{t)\^ + 2uj \\a{T)\\jjrdT <\a{to)\^ + 2 f {f,a{T))dT, 

J to J to 

for all <tQ<t. 

We write a solution a{t) to (13.11 ) in the form 

(3.3) aj{t) =aj + bj{t), 

where a is a fixed point whose properties were exhibited in Section 2. We now 
show that this fixed point is the exponential global attractor. In particular, the fixed 
point is unique. 

Theorem 3.4. Let a £ be a fixed point of (13.11) for c G (3/2, 5/2] and a{t) be a 
solution with a(0) G andaj{0) > Ofor all j. Then 

(3.4) !6(t)|2 < 15(0)126-27^^*. 
Proof. As before we write 

= 2-^■'■/='2^/Vo^^^i• 
Now let Bj{t) := 2'=j/36j(t). Then 

(3.5) aj{t) = aj + 2-''i'^Bj{t). 



Then the system (13.11) reduces to 
(3.6) 



g 
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where j > and A^i = B^i = 0. Following the procedures in the inviscid case 

given in Q we multiply (13.61) by Bj and sum to obtain 

(3.7) 

Since Ak+Bk{t) > for all i > and lim,fc^oo ^fc = 0, we have that lim inf yfc_»oo Bk+i{t) > 
for all t > 0. Then due to the fact that Aj decreases super-exponentially (see 
Lemma |2?T]) . and Yl'jLo'^^'' -^j integrable, we can use the dominated conver- 
gence theorem to obtain 

oo oo 

(3.8) fo'^''2'/^ Yl ^'^"'^^B.itf - /o"'/'2^/6 ^ 2-2-^/35,(0)2 
i=o j=Q 



/■ r 

where j > and A_i = = 0. Then we have 
(3.9) 

oo 
j=0 

oo oo 



dr, 



oo 



j=0 



j=Q 
oo 



1/2 



j=0 

oo 



+ ^j+l^j 



1/2 .1/2 ^ ^3^J+l 



■j+2 



j=0 

where we used equation (12.11) in the last inequality. Now using Cauchy-Schwarz 
inequality and Lemma [241 we obtain 

(3.10) 

oo oo 

E H^-i^.-i^. - AjB.+iB, - A,^,B] - ^2^^i?|] < E B]. 

j=0 j=0 

Therefore 
(3.11) 

OO OO oo 

E 2-2-^/3b, (t)2 - ^ 2~^^^/'B,{0f < 2/o^/22~c/6^^ / ^ 2«i?, (r)2 dr 

„■ n n n 
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Hence for bj = l-^^^^Bj we have 



j=0 j=0 j=0 

(3.12) 

E 

j=0 



which imphes that 

(3.13) \b{t)f < |6(0)pe-2^^*. 

4. Dissipation Anomaly 



□ 



Here we study the energy dissipation in the limit of vanishing viscosity. For 
convenience, solutions to the dyadic model with viscosity > will be denoted 
by a'^{t) in this section. The fixed point (which is unique in both viscous and 
inviscid cases) will be denoted by a^. Now given a solution a^{t) to the inviscid 
dyadic model, we define its anomalous energy dissipation rate as follows 

e,o(t) := (a0(t),/)-i||a0(t)|2 

in the sense of distributions. Due to the energy inequality, e^o > and hence e^o 
is a Borel measure. The following theorem was proved in [?]: 

Theorem 4.1. Let a^{t) be a solution to the inviscid dyadic model on [0, 00). Then 

rT 



1 

- deAt) =■■ ed > 0. 



lim 

T^oo T 

Since the global attractor of the viscous model A"^ converges to the global at- 
tractor of the inviscid model as v ^ (see Lemma 123] ) . we have the following 

Theorem 4.2. Let a^{t) be a solution to the viscous dyadic model on [0, 00). Then 



lim i / u\\a''{t)\\]j,dt^eA>Q, 



as u ^ 

Proof. Due to the energy inequality we have 

1 1 1 f*+^ 1 r*+^ 

l^V{t)?-,^y{t+T)\^ < -u- \\a^^s)\\l,ds+- {a^s),f)ds. 

Hence, 

limsupi r u\\a''{t)\\l^dt< lim i {a'' (s), f) ds = {a-^ , f). 

On the other hand, note that the fixed point a*^ (since it is a regular solution) satis- 
fies the energy equality 

- [a ,/). 
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Now for any 6 > 0, there exists A^, such that 

N 
3=0 

Since a'^{t) in t^, we have 



hminf ^ / v^i^^d'itf^dt > > - ^■ 

Therefore, 

liminf;^ / u\\a''{t)\\jj,dt>u\\a''\\j^i = {a'',f). 
Since {a'^, f) {a^, f ) as u 0, we obtain 

lim 1 f y^ml, dt = (a",/) = a[]/o =: ea > 0. 



□ 



5. Dissipation length scale 

Kohnogorov's theory of turbulence predicts that the energy density in the inertial 
range is 

(5.1) £{k) ^ e'J\-'/^ 

followed by a rapid decay after the dissipation wave number 



1 

Cd \ 4 



(5.2) Kd 

For the dyadic model the energy density is defined as £'(2-^) = a^2~-'. In the 
inviscid case one can easily check that the energy density for the fixed point (which 
is a global attractor) is 

(5.3) £iK)^efK-^'/-'-\ 

By Lemma [23] the energy density on the global attractor of the dyadic system with 
small positive viscosity is close to ( 15.31 ) in the inertial range. Moreover, Lemma lZT] 
can be used to determine the dissipation wavenumber for the model 



(5.4) Kd := 2 



J 



1 2 

^3/2 \ 4 ■ 3-c 1 2 



Lemma 2.1 is valid for all c in our range of interest, including c = 1. The inequal- 
ity (2.2) implies that ultimately the energy in each shell decays very rapidly with 
increasing j. When we invoke monotonicity of the sequence Aj, (2.2) ensures that 
this very rapid decay occurs at the dissipation wave number scale given by (5.4). 
However for technical reasons we require c > 3/2 in order to prove monotonicity. 
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As we discussed in Section 1, the appropriate range for c, that arises from Bern- 
stein's inequality applied to estimate the nonlinear term in the Navier-Stokes equa- 
tions, is c E [1, 5/2]. We note in the end point case c = 1 the expressions for the 
energy density and the dissipation wavenumber (15.31 ) and (15.41 ) coincide with (15.11) 
and (1521) . i.e., 

8{K)^eT^-"\ for ^=1- 

At the end point value c = 5/2, which corresponds to complete saturation of Bern- 
stein's (and Sobolev) inequalities, we have 

~ e^/\-8/3^ for c = 5/2. 
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